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1 .  Introduction 

Imaging  of  single  atoms  in  the  field  ion  microscope  is  believed  to  occur 
as  a  result  of  local  field  enhancement  around  kink  sites  and  around  single 
metal  atoms  on  flat  crystal  planes.  These  local  electric  fields  of  the 
order  of  V/X  arise,  via  Poisson's  equation,  from  local  enhancements  of  the 
surface  electron  density.  To  understand  this  phenomenon  from  a  microscopic 
point  of  view,  it  is  mandatory  to  calculate  the  electronic  charge  distribu¬ 
tion  and  the  resulting  electrostatic  field  self -consistently. 

Classical  electrostatic  theory  predicts  that  the  electrostatic  potential 
around  a  hemispherical  protrusion  of  radius  R  on  an  otherwise  flat  and 
structureless  metal  surface  is  given  byC1) 


<t>  =  —4xa0z  (1  -  R3/r3 )  (1) 

where  ct0  is  the  excess  charge  far  from  the  protrusion  and  r  is  measured  from 
the  center  of  the  protrusion.  locally  the  excess  charge  redistributes 
itself  into 

a  -  a0(l  -  R3/r»)  (2) 

on  the  plane,  and 

a  =  3<t0z/R  (3) 

on  the  sphere.  We  note  that,  at  the  apex  of  the  boss,  the  electric  field  is 
three  times  its  value  at  infinity. 
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The  first  quantum  mechanical  calculation  of  the  self-consistent  charge 
density  and  electric  field  distribution  at  a  metal  surface  was  performed  by 
Lang  and  Kohn(2)  for  a  structureless  jellium  surface  using  the  local-den¬ 
sity  approximation  to  density  functional  theory.  A  systematic  study  of 
this  problem  was  recently  performed  by  Gies  and  Gerhardts  (2 )  and  also  by 
Schreier  and  Rebentrost (4) .  Including  the  crystal  structure,  Ihglesfield 
has  performed  self-consistent  field  calculations  for  A1  and  Ag  sur¬ 
faces  (5'6).  (For  a  recent  review  of  field  adsorption,  field  desorption  and 
field  evaporation,  see  ref .7.) 

To  understand  field  effects  in  the  field  ion  microscope,  one  would  like 
to  calculate  the  local  field  at  kink  sites  or  around  single  atoms  on 
densely  packed  planes  of  transition  metals.  As  a  model  of  the  latter,  we 
will  report,  in  this  paper,  self-consistent  calculations  of  the  electric 
field  around  a  metal  atom  adsorbed  on  a  structureless  jellium  surface  using 
the  local  density  functional  theory.  As  examples,  we  consider  a  light  and  a 
medium  heavy  metal  atom,  Ti  and  Nb,  respectively,  on  a  jellium  metal  with  a 
Wigner-Seitz  radius  rs=3.0  bohr.  We  will  restrict  ourselves  to  the  polar¬ 
ity  where  the  field  points  away  from  the  surface,  appropriate  for  the  field 
ion  microscope.  Although  our  model  describes  the  chemisorption  of  metal 
atoms  on  a  free  electron  metal  in  the  presence  of  a  field,  we  will  see  below 
that  thQ  results  have  quite  general  features  that  allow  a  discussion  of 
field  effects  in  chemisorption  on  transition  metals. 

In  the  next  section  we  collect  the  relevant  equations  and  specify  the 
boundary  conditions.  In  vbe  following  sections  we  discuss  numerical  exam¬ 
ples,  i.e.  the  adsorption  of  Nb  and  Ti  atoms  on  jellium  in  the  presence  of  an 
external  field.  We  will  recover  the  field  enhancement,  described  classi¬ 
cally  by  (1-3),  albeit  of  smaller  magnitude.  We  will  also  produce  the 
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adsorption  potential  well  for  these  atoms  and  calculate  the  evaporation 
field  strength  above  which  surface  atoms  are  no  longer  bound  to  the  sur¬ 
face.  Despite  the  simplicity  of  the  jellium  model  we  find  very  good  agree¬ 
ment  with  experimental  values,  and  also  with  earlier  semi-empirical  calcu¬ 
lations  based  on  the  ASED-MO  method (7/8).  We  will  also  discuss  the  dipole 
moment  and  the  polarizability  of  the  adsorbed  atom  to  make  a  connection 
with  some  experimental  results. 

2 .  The  Model 

In  the  introduction  we  have  described  electric  fields  by  using  classical 
electromagnetic  theory,  which  assumes  that  the  surface  of  a  metal  is  a 
mathematical  plane  with  excess  charges  and  a  dipole  layer  at  which  the 
normal  component  of  the  electric  field  drops  discontinuously  to  zero,  at 
least  for  a  perfect  conductor.  On  real  surfaces,  however,  the  electron 
distribution  and  electric  fields  vary  smoothly  over  distances  of  a  few  ang¬ 
stroms.  A  simple  model (2/9)  that  exhibits  these  features  is  the  jellium 
model  of  a  metal,  in  which  it  is  assumed  that  the  ionic  lattice  can  be  smo¬ 
othed  into  a  uniform  positive  charge  density  n+  that  drops  to  zero  abruptly 
half  a  lattice  constant  above  the  topmost  layer  of  ion  cores.  In  front  of 
this  jellium  metal  we  add  a  metal  atom(10)  around  which  we  will  study  the 
local  electron  and  field  distribution  within  the  framework  of  density  func¬ 
tional  theory.  The  latter  determines  the  exact  ground-state  electron 
(number)  density  as  the  self-consistent  solution  of  the  equations 

n(r)  =  -5-  Y  6  )  |  i>i/  (*j|  2 

v 


(4) 
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2^-V2  +  veff  [n;r]  -eu 


4>v  (r)  =  0 


(5) 


In  the  local-density  approximation,  the  effective  potential 


veff  [n'r]  =  <t>  (r)  +  Mxc  (n(r))  (6) 

is  the  sum  of  the  electrostatic  potential  as  seen  by  an  electron,  deter¬ 
mined  from  Poisson' s  equation 


v*<t>  (r)  =  4jre[n+  -  n(r)]  (7) 

and  the  exchange  and  correlation  contribution,  nxc ,  as  discussed  in  refs.  2 
and  9.  e>0  is  the  elementary  charge. 

To  study  chemisorption  of  an  atom  on  jellium,  one  rewrites  (5)  in  terms  of 
a  Lippmann-Schwinger  equation  (10) 

^,MA  (r)  =  (r)  +  Jd3^  gM  (r^  )5veff  ( r'  )^MA  (r' )  (8) 

where  is  the  solution  of  (5)  in  the  absence  of  the  adsorbed  atom,  is 
the  Green' s  function  of  the  bare  metal  and  6veff  is  the  difference  of  the 
effective  potential  (6)  with  and  without  the  adsorbed  atom.  Because 
6vef f  (r)  approaches  zero  rapidly  outside  the  adatom,  one  can  restrict  the 
integral  in  (8)  to  a  sphere  centered  at  the  adatom,  bounded  by  a  surface  S 
outside  of  which  5veff  is  negligible. 

To  include  an  external  electric  field,  one  must  add  an  excess  surface 
charge  to  the  metal.  This  is  done  first  in  the  bare-metal  problem  by  impos- 
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ing  the  condition 


a0  =  e 


OO 


(n+9(-z)  -  nM  (z,F0))dz  =  F0/4tt 


-OO 


(9) 


on  the  charge  density  of  the  bare  metal,  leading  to  wavefunctions  (r,F0) 
and  electron  densities  rF1  (r,F0)  that  depend  explicitly  on  the  field 
strength,  F0,  far  from  the  surface.  In  the  presence  of  the  adatom  one  then 
requires  that  outside  the  sphere  S 

nMA  (r;  F0)  -  nM  (r;  F0)  for  |  r|  -  oo  (10) 

which  acts  as  a  boundary  condition  on  the  solution  of  (8).  Thus,  the  effect 
of  the  external  electric  field  is  imposed  in  the  metal  surface-adatom 
system  through  the  condition  (9). 

While  the  external  electric  field  does  not  change  the  structure  of  the 
one-electron  equations,  it  does  affect  the  numerical  procedure.  In  the 
absence  of  external  electric  fields,  the  effective  potential  of  a  bare 
metal,  veff  [nM  ;  z] ,  will  approach  a  constant  for  z-*oo,  usually  taken  as  the 
zero  of  energy.  However,  in  the  presence  of  an  external  electric  field  of 
asymptotic  strength  F0,  veff[nM;z]  will  grow  as  eF0z  for  large  z.  It  is 
then  advantageous  to  choose  the  Fermi  energy  as  the  energy  reference,  ep  = 
0.  With  this  choice  the  energies  of  all  occupied  levels  are  negative. 

To  solve  (8),  one  represents  the  wavefunctions  with  indices  E,  m  and  1, 
t/)£ml  i  as  r-dependent  linear  combinations  of  free-particle  solutions,  fol¬ 
lowing  ref.  10, 
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Eml  (r> 


V 

r  =o 


cEmir  (r)3l'  (Pr)  +  sEmir 


(r)h+ 


>] 


(pr)  Yrm(Q) 


(11) 


Here  p  =  %/eT and  and  h^,  are  the  spherical  Bessel  functions.  A  complica¬ 
tion  arising  in  the  presence  of  an  electric  field  is  the  fact  that  the 
potential -energy  barrier  of  field  adsorption  varies  from  several  eV,  i.e. 
the  cohesive  energy  in  the  absence  of  the  field,  down  to  zero  at  the  evapo¬ 
ration  field  strength.  Recalling  that  the  theory  calculates  the  energy 
difference  AE^-  between  the  bare  metal  surface  and  the  metal  surface-adatom 
system,  i.e.  the  total  energy  of  the  adatom  of  the  order  of  105  eV  for  Ti  and 
Nb,  we  need  a  relative  precision  for  the  total  energy  difference  of  the 
order  of  about  10-7 . 


3 .  Binding  Energies  and  Energy  Barriers 

We  now  present  numerical  results  for  two  systems,  titanium  and  niobium, 
respectively,  on  a  jellium  with  rs=3.0  bohr. 

In  Fig.l,  we  show  the  adiabatic  energy  curves  for  niobium  on  jellium  for 
several  field  strengths.  In  zero  field  the  depth  of  the  surface  potential 
is  5.2eV,  which  is,  not  surprisingly,  somewhat  less  than  the  cohesive  energy 
of  Nb,  7.47eV,  because  the  jellium  lacks  p-  and  d-orbitals.  This  deficiency 
is  also  the  cause  of  the  rather  weak  repulsion  at  short  distances.  As  we 
apply  a  field,  we  note  a  shift  of  the  potential  minimum  by  about  0.2  X  away 
from  the  surface.  We  will  trace  this  shift  to  a  transfer  of  electronic 
charge  from  outside  the  adatom  to  the  region  between  the  adatom  and  the 
surface.  With  a  field  applied,  the  ground  state  energy  curves  must  assume 
an  asymptotic  form  -eF0z  for  large  z,  appropriate  for  a  singly  charged  posi¬ 
tive  ion.  Ionization  of  the  adatom  can  occur  when  its  highest  occupied 
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level  is  lifted  by  the  field  energy  term,  eFcz,  above  the  Fermi  level  of  the 
jellium.  If  the  ionization  level  were  not  shifted  and  broadened  by  the  in¬ 
teraction  with  the  metal,  the  changeover  to  the  asymptotic  form  would 
happen  abruptly  at  the  apex  of  the  potential  energy  curve,  i.e.  at  the  point 
where  the  diabatic  energy  curves  for  the  neutral  and  ionic  spedies  cross. 
For  the  fields  chosen  for  Fig.l,  the  apex  is  so  close  to  the  metal  surface 
that  considerable  interaction  between  the  adatom  and  the  metal  is  still  in 
effect.  This  results  in  a  considerable  broadening  of  the  ionization  level 
of  the  adatom,  so  that  only  partial  charge  draining  occurs  in  the  apex 
region . 

We  will  call  the  energy  difference  between  the  minimum  of  the  surface 
potential  and  the  local  maximum  at  its  apex,  the  activation  energy,  Q(F0), 
for  ionization.  In  zero  field,  the  activation  energy  is  equal  to  the  binding 
energy  of  the  adatom.  We  note  that  the  activation  energy  becomes  zero  for 
the  evaporation  field  strength  which  we  estimate  for  Nb  on  jellium  to  be 
3.6V/A,  which  compares  very  favorably  with  the  experimental  value  for  Nb, 
Fev=3*5V/iL  Similar  results  are  obtained  for  Ti.  Binding  and  activation 
energies  are  collected  in  Table  1. 

It  has  been  suggested  (8)  that  the  field  dependence  of  the  activation 
energy  for  different  metals  obeys  a  universal  scaling  law  if  one  plots 
Q(F0)/Q(F0=0)  as  a  function  of  f=F0/Fev.  In  a  simple  model,  this  scaling  law 
is  given  by 

Q(Fo)/Q(F0=0)  -  VTT  +  <n[(l-^H)/(l+^hf)]  (12) 
Experimental  data  on  tungsten C11)  and  theoretical  results  obtained  by  the 
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ASED-MO  method (8)  have  confirmed  this  conjecture,  as  do  the  present  results 
for  Ti  and  Nb. 

4.  Charge  Density  and  External  Electrostatic  Field 
In  Fig. 2  we  show  the  electron  density  shift 

*n(r,  F0 )  =  n(r,  F0)  -  n(r,  F0=0)  (13) 

due  to  the  applied  field,  demonstrating  an  overall  charge  transfer  from  the 
region  outside  the  adsorbed  atom  into  the  region  between  the  atom  and  the 
bulk  and  further  into  the  bulk.  From  (13)  we  can  calculate,  via  Poisson's 
equation  (7),  the  external  electrostatic  potential  as  the  difference  of  the 
total  electrostatic  potential  in  the  presence  of  an  excess  charge  of  den¬ 
sity  <70=F0/4rr ,  and  in  the  absence  of  it, 

^ext(r)  =Ves(r'F0)  -  Ves(r,F0=0)  .  (14) 

Equipotential  contours  are  plotted  in  Fig. 3  for  Ti  and  Nb.  The  expulsion  of 
the  external  field  from  the  region  of  the  metallic  adatom  is  beautifully 
demonstrated,  the  effect  being  larger  for  the  heavier  Nb  atom.  To  compare 
this  quantum  mechanical  result  with  the  classical  theory  contained  in 
(1-3),  we  have  approximated  the  adsorbed  atom  by  a  hemispherical  boss  with 
the  radius  chosen  in  such  a  way  that  the  zero  equipotential  lines  coincide 
at  the  apex.  The  result  of  this  classical  fit  is  also  shown  in  Fig. 3. 
Except  on  the  sides  of  the  adatoms,  the  agreement  between  the  quantum 
mechanical  and  classical  theories  is  quite  good,  at  least  at  the  semi- quan¬ 
titative  level.  There  are  some  differences,  of  course,  if  we  look  at  deta- 
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ils,  which  we  do  in  Fig. 4,  where  we  plot  the  external  field  along  a  line 
through  the  apex  of  the  adatom.  Compared  with  the  field  in  the  absence  of 
the  adatom,  we  again  note  the  expulsion  of  the  field  from  the  adatom  region 
which  results  in  an  enhancement  of  the  field  just  outside  the  adatom.  We 
note,  however,  that  this  field  enhancement  at  the  apex  is  not  as  much  as 
classical  theory  predicts,  i.e.  a  factor  3.  Rather,  the  partial  field  pene¬ 
tration  into  the  adatom  results  in  a  smearing  out  of  the  field  as  a  reflec¬ 
tion  of  the  adjustability  of  the  electronic  distribution  at  the  surface, 
reducing  the  enhancement  effect.  I.ote  again,  that  for  the  heavier  adatom, 
Nb,  one  is  closer  to  the  classical  result.  It  is  also  noteworthy  that  the 
Friedel  oscillations  inside  the  bulk  jellium  are  quite  similar  with  and 
without  the  adatom. 

5 .  Core  states 

As  we  saw  in  Fig. 4,  there  is  a  substantial  expulsion  of  the  electrostatic 
field  from  the  region  of  the  adatom.  To  understand  the  effect  of  an  exter¬ 
nal  electrostatic  field  on  the  adatom,  we  start  from  a  bare  metal  surface 
with  an  external  potential  (r;F0)  due  to  a  surface  charge  o0  produc¬ 
ing  an  asymptotic  field  F0  ^47r<70 .  When  we  bring  an  atom  into  this  field,  it 
will  be  polarized,  producing  an  induced  field  with  a  potential  (r)  which 
cancels  the  original  field  to  a  large  extent  in  the  vicinity  of  the  atom. 
Because  the  adatom  is  in  close  contact  with  the  metal  surface,  there  will, 
in  addition,  be  a  charge  transfer  from  the  atom  to  the  surface  involving 
valence  electrons  of  the  adatom.  This  charge  transfer  also  affects  the 
core-level  states,  as  demonstrated  in  Fig. 5  where  we  plot  the  energy  eigen¬ 
value  of  the  Is  level  of  Ti  as  a  function  of  the  nuclear  distance  from  the 
jellium  edge  for  different  field  strengths.  As  discussed  in  ref. 10,  the 
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energy  eigenvalues,  calculated  within  the  framework  of  density  functional 
theory,  do  not  directly  give  the  energy  to  remove  an  electron  from  the 
respective  states,  however,  their  distance  dependence  reflects  the  varia¬ 
tion  of  the  chemical  shift.  Close  to  the  surface  there  is  a  net  transfer  of 
electrons  from  the  metal  to  the  adatom  generating  more  repulsion  which  in 
turn  causes  the  Is  level  of  the  adatom  to  be  higher  than  that  of  the  iso¬ 
lated  Ti  atom.  As  we  move  the  adatom  away  from  the  surface  there  is  even¬ 
tually  some  charge  draining  to  the  surface  creating  more  attraction  so  that 
the  Is  level  decreases.  Even  at  4.0  bohr  it  is  still  above  the  Is  level  of 
the  Ti+  ion. 

This  effect  suggests  a  way  to  assign  a  partial  charge,  q,  to  the  adatom 
without  the  inherent  arbitrariness  connected  with  the  Mulliken  population 
analysis.  In  lowest  order,  we  set 

q/e  =  [E (Is)  -  E0(ls)]/[E+(ls)  -  E0(ls)j  (15) 

where  E(ls),  E0  (Is)  and  E-|_(ls)  are  the  energies  of  the  Is  level  of  the 
adatom,  the  isolated  atom  and  the  isolated  ion,  respectively.  This  quantity 
is  also  plotted  in  Fig. 5  as  a  function  of  the  adatom  position.  We  note  that 
at  4.0  bohr  these  partial  charges  agree  with  the  slopes,  -qF0 ,  of  the 
ground-state  energy  curves. 

6.  Dipole  Moment 

Experiments  have  been  devised  to  extract  the  dipole  moment  and  the  pol¬ 
ar  is  abil  ity  of  a  metallic  adatom  from  field  emission  work  function 
changes  (12) .  The  former  quantity  can  be  rigorously  defined  in  our  theory  as 
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M(F0)  =  -e/d3r  r[n(MA)  (r,F0)  -  n(M)  (r,F0)  -  Z6(r-r0)]  (16) 

where  n^A)  and  n(M)  are  the  electron  number  densities  with  and  without  the 
adatom,  respectively,  and  Z  is  the  nuclear  charge  of  the  adatom  at  position 
r0.  We  note  that  the  total  charge  in  the  two  situations  is  identical,  i.e. 
producing  an  asymptotic  field  F0 .  This  guarantees  that  the  definition  (15) 
of  the  dipole  moment  does  not  depend  on  the  origin  chosen. 

Looking  next  at  the  dependence  of  the  dipole  moment,  n,  of  adsorbed  Ti  on 
the  adatom  position,  d,  one  would  expect  a  linear  rise  of  n  with  d,  if  the 
charge  transfer  from  the  adatom  to  the  metal  were  independent  of  the 
adatom  position.  However,  one  finds  that  the  rise  is  less  than  linear  due 
to  partial  re-arrangement  of  charges  as  the  adatom  is  pulled  away  from  the 
surface.  Obviously,  once  one  is  far  enough  beyond  the  hump  of  the  Schottky 
barrier,  the  adatom  gets  ionized  and  a  linear  rise  in  ^  should  persist  up  to 
the  point  where  double  ionization  becomes  possible. 

In  the  experimental  literature  one  parametrizes  the  dipole  moment  as 

M(F0)  =  /x(F0=0)  +  qF0  (17) 

in  terms  of  the  asymptotic  field  F0,  defining  the  (field-dependent)  polari- 
sability  of  the  surface  atom  q(F0),  rather  than  introducing  hyperpolariza¬ 
bilities,  gradient  terms  etc.  Values  of  /x(F0)  and  a(F0)  are  given  in  Table  2. 
Measured  polarisabilities(12)  at  field  strengths  F0<0.4V/X  decrease  mono- 
tonically  from  11.94&3  for  Ta  to  2.74&3  for  Pt.  Our  numbers  for  Ti  and  Nb 
are,  indeed,  close  to  that  for  Ta. 

The  field-induced  dipole  moment  of  a  metal  atom  at  a  kink  site  or  on  top 
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of  a  closed-packed  plane  has  been  used  by  Tsong  and  Muller  (13)  to  estimate 
the  local  field  enhancement 


F(r)  =  F0  +  a(F0) 


3n(F0-n)  -  F0 

r3 


(18) 


Here  F0  is  the  constant  electric  field  far  from  the  surface  and  n  is  a  unit 
vector  in  the  direction  of  r.  Along  a  line  perpendicular  to  the  jellium 
surface  and  through  the  center  of  the  adatom,  the  field  from  (18)  is  given 
by 


FZ  (z)  =  F0 


1  + 


2a  (F0) 
z3 


(19) 


whereas  classical  theory  gives,  from  (1-3), 

Fz  (z)  =  F0 

Tsong  and  Muller(13)  measured  r  in  (18)  from  the  nucleus  of  the  adatom, 
rather  than  from  the  center  of  mass  of  the  charge  distribution.  This  over¬ 
estimates  the  field  enhancement  considerably.  We  will  present  a  detailed 
discussion  of  the  effects  of  field  enhancement  on  field  adsorption  else¬ 
where  . 


1  + 


2R3 


(20) 


7.  Conclusions 

We  have  presented  self-consistent  calculations  of  the  electronic  structure 
of  Ti  and  Nb  atoms  adsorbed  on  a  jellium  metal  in  the  presence  of  an  exter- 
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nal  electric  field.  We  have  calculated  ground-state  energy  curves, 
extracting  activation  energies  for  field  evaporation  and  evaporation  field 
strengths,  finding  good  agreement  with  data  for  transition  metals.  Like¬ 
wise,  the  calculated  dipole  moment  and  polarisability  of  these  metallic 
adatoms  are  consistent  with  experimental  values  on  other  metals. 

The  field  enhancement  at  kink  sites,  along  steps,  and  above  atoms  in  open 
surfaces,  plays  a  crucial  role  in  the  image  formation  in  the  field  ion 
microscope,  as  demonstrated  by  the  fact  that  densely  packed  surfaces 
cannot  be  resolved.  Our  results  show  that  a  field  enhancement  by  a  factor 
1.5  to  2  must  be  expected  1-2&  above  an  isolated  metal  atom.  The  implica¬ 
tions  of  this  effect  for  field  adsorption  of  rare  gases  will  be  explored  in 
a  forthcoming  paper. 
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Table  1:  Energy  barriers  Q(F0)  (in  eV) .  The  zero  field  value,  Q(0),  is  the 
binding  energy. 


adatom 

F(V/A) 

0.0 

2.2 

2.6 

3.0 

3.4 

Ti 

4.3 

0.57 

0.25 

0.09 

- 

Nb 

5.2 

1.3 

0.64 

0.20 

0.05 

Experimental  results (14) : 


binding  energy 
evaporation  field 


4 .85  eV  for  Ti  and  7 .47  eV  for  Nb 
2.5  V/%.  for  Ti  and  3.5  V/%  for  Nb 
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Table  2:  Dipole  moment,  n  (in  Debye)  and  polarizability,  a  (in  %?),  at  the 

equilibrium  distance,  as  a  function  of  field  strength. 


adatom 

F(V/&) 

2.2 

2.6 

3.0 

3.4 

Ti 

A* 

7.3 

10.3 

12.5 

- 

a 

9.9 

11.9 

12.4 

- 

Nb 

7.3 

11.1 

13.1 

15.0 

a 

10.1 

12.8 

13.1 

13.3 
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Figure  Captions 

Fig.l:  Potential  energy  for  Nb  on  jellium  with  rs=3.0  bohr  as  a  func¬ 

tion  of  nuclear  distance,  z,  from  the  jellium  edge.  F0  in  V/iL 

Fig. 2:  Electronic  density  shift  in  an  electric  field  for  the  adatom  at 

the  minimum,  d,  of  the  potential  energy.  Left:  Ti  (rs=3.0  bohr 
at  d=2.5  bohr)  in  F0=2.6V/il.  Right:  Nb  (rs=3.0  bohr  at  d=2.75 
bohr)  in  F0=3.0V/JL  Contours  Sn  =  0,  ±0.0001,  ±0.002,  ±0.005, 
±0.001,  ±0.002,  ±0.005  a.u.  Solid  (broken)  lines  non-negative 
(negative)  values.  Contours  in  the  core  region  are  omitted. 

Fig. 3:  External  electrostatic  potential,  upper  panels  for  the  same 

systems  as  Fig. 2;  lower  panels:  classical  results.  Contour 
lines  -0.05,  -0.02,  -0.01,  0.0,  0.03,  0.1,  0.2,  0.3,  0.4,  0.5,  0.6 
Ry.  Solid  (broken)  lines  non-negative  (negative)  values. 

Fig. 4:  Electrostatic  field  strength  along  line  through  the  center  of 

the  adatom  for  the  same  systems  as  in  Fig. 2.  Curves  A  and  B: 
without  and  with  the  adatom;  curve  C:  classical  result. 

Fig. 5:  Is  energy  levels  and  estimated  charge,  q,  for  Ti.  The  levels  in 

the  isolated  atom  and  ion  are  indicated  on  the  right  of  the  left 
panel . 
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